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ABSTRACT 

We study analytically a gravitational lens due to a deformed star, which is modeled by 
using a monopole and a quadrupole moment. Positions of the images are discussed for 
a source on the principal axis. We present explicit expressions for the lens equation for 
this gravitational lens as a single real tenth-order algebraic equation. Furthermore, we 
compute an expression for the caustics as a discriminant for the polynomial. Another 
simple parametric representation of the caustics is also presented in a more tractable 
form. A simple expression for the critical curves is obtained to clarify a topological fea- 
ture of the critical curves; the curves are simply connected if and only if the distortion 
is sufficiently large. 
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1 INTRODUCTION 

Some recent highly accurate numerical simulations of ro- 
tating compact objects such as a neutron star elucidate 
many features of compact objects such as an equilibrium 
configuration and its stability (e.g., Shibata et al. 2002, Ster- 
gioulas et al. 2004, Zanotti and Rezzolla 2002, Zdunik et al. 
2004). The deformation of their configurations causes a non- 
spherical distortion of the surrounding spacetime, which has 
not been observed yet. Probing this warp of the spacetime 
provides a strong verification of the theory of the general 
relativity (Will 1993) and a deep insight into the relevant 
physics especially about the nuclear matter (e.g., Shapiro 
and Teukolsky 1983), because the deformation depends not 
only on the gravity but also on the equation of state of the 
internal matter. 

One possibility is to study effects of the deformation 
on the light propagation near the compact object. Before 
doing numerical studies in detail, it is important to find 
some physical properties by using a simple analytic model. 
One reason is that one needs many light ray tracings in 
numerical investigations, and consequently it is difficult to 
clarify quantitatively the parameter dependence. 

We perform the multipole expansion of the external 
gravitational field, which can be assumed to be weak along 
the line of sight even if the internal gravity is strong. Namely, 
one can consider only the linearized potential. This does not 
mean that our system is linear, because the lens equation is 



* Email: asada@phys.hirosaki-u.ac.jp 



non-linearly coupled with respect to the image positions as 
shown in the next section. 

The quadrupole moment depends on the rotation speed 
and the equation of the state for the internal matter. Many 
years ago some comprehensive studies had been done for 
equilibrium configurations of rotating stars in the Newtonian 
gravity, some of which are summarized in a famous textbook 
(Chandrasekhar 1987). Numerical investigations on this is- 
sue in the general relativity are still going on as mentioned 
in the beginning of this section (e.g., Shibata et al. 2002, 
Stergioulas et al. 2004, Zanotti and Rezzolla 2002, Zdunik 
et al. 2004). In an extreme case called as rapid rotation, 
the quadrupole deviation of order of 10 -3 may be possible 
even for a neutron star. These results show that the defor- 
mation is dominated by the quadrupole moment and higher 
moments such as the octapole moment are much smaller. 
In addition, a contribution of the octapole moment in the 
external field decays faster than that of the quadrupole mo- 
ment, because the ratio between the octapole potential and 
the quadrupole potential is of the order of the ratio of the 
stellar radius to the distance between the star and the field 
point, which becomes the impact parameter in the gravi- 
tational lensing. This ratio is so small that one can ignore 
approximately the octapole and higher moment in the lens- 
ing study. Furthermore, radio observations such as the VLBI 
will be improved. If the accuracy is improved sufficiently, the 
present result will have a direct relevance with such an ob- 
servation. Hence, we consider in this paper the truncated 
multipole expansion up to the quadrupole moment. 

This consideration in the gravitational lens is not new 
but rather classical though it was done in a different context. 
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For instance, the light deflection caused by the quadrupole 
moment of the Sun was discussed in the context of tests of 
the general relativity (Epstein and Shapiro 1980). Except 
for a case of the source on the equatorial plane, we need 
solve the lens equation numerically. A spherically symmetric 
galaxy and a quadrupole tidal part were extensively studied 
by Kovner (1987), in which most models are an extended 
sphere and the quadrupole part, and another model is a 
star with quadrupole shear induced by a galaxy where the 
star is located. The last model seems similar to our system 
of a deformed star. However, it should be noted that this 
extra shear does not come from a potential of the star but 
from that of the galaxy. In other words, the dependence of 
the shear on the radial coordinate is different from that for 
the deformed star. 

It is believed that gravitational lensing study needs 
some numerical techniques except for a few special cases 
such as a point mass lens or an axisymmetric mass distri- 
bution along the line of sight, because the lens equation is 
nonlinearly coupled (for instance, Schneider et al. 1992). In 
recent, however, some development has been achieved: It 
has been shown that for two cases of a binary lens (Asada 
2002) and an ellipsoidal lens (Asada et al. 2003) the lens 
equation is reduced to a single polynomial. An advantage of 
this approach is saving time and computer resources, and 
another one is to enable us to study analytically the lensing, 
for instance to find out an analytic expression for caustics 
(Asada et al. 2002, 2003) which are curves on the source 
plane, and the number of images changes if a source crosses 
the caustics. This expression for the caustics is nothing but a 
discriminant for the polynomial, which determines a change 
in the number of roots for the equation: As known in algebra, 
the discriminant vanishes in a case of multiple roots, which 
corresponds to merging images or giant arcs in the gravita- 
tional lens. Classification of the critical curves and caustics 
is of great importance. For instance, it was done intensively 
for a binary gravitational lens (e.g., Dominik 1999). 

Furthermore, it has been shown also that one can de- 
rive a single polynomial equation for any set of the lens 
equations in polynomials by using the Euclidean algorithm 
(Asada et al. 2004). In practice, however, in a case of the 
lens equation which is higher than fifth degree, this alge- 
braic procedure needs a computer with a fast CPU and a 
huge memory, and furthermore the resultant equation may 
become too complicated to handle any more in the Carte- 
sian coordinates. In order to avoid this problem, the polar 
coordinates are adopted in this paper so that a simple result 
can be obtained. 

This paper is organised as follows. First, we show that 
the lens equation is reduced to a single polynomial equation. 
Next, the analytic expression for the caustics is computed 
as the discriminant for the polynomial. Finally, a paramet- 
ric representation for the critical curves and caustics are ob- 
tained to clarify a topological feature of the critical curves. 



2 LENS EQUATION FOR A LENS WITH 
QUADRUPOLE MOMENT 

2.1 Lens Equation 

The lens equation relates the image angular position vector 
to the source angular position vector /3 as 



/3 = 



D 



LS 



(1) 



where we used the thin-lens approximation, Dg is the an- 
gular distance from the observer to the source, and £>Lg i g 
that from the lens to the source. The vector £, which de- 
notes the impact parameter to the light ray, is related to 
the image position as £ = ^L^' wnere ^ s ^ ne an g mar 
distance from the observer to the lens. 

Up to quadrupole moments in multipole expansions, the 
deflection angle a — (a 1 , a 2 ) is given by 



j,k=l j=l 



Qij 



•(2) 



Here M denotes the mass of the lens, and the quadrupole 
moment tracefree on the lens plane is defined as 



Qi. 



p(XiXj 



X\ 2 )d 3 X, 



(3) 



where p denotes the mass density. It is worthwhile to men- 
tion that there are no contributions from quadrupole mo- 
ment whose components are along the line of sight (Asada 
and Kasai 2000). 

Now, a renormalisation is done in units of the Einstein 
ring angular radius defined as (for instance, Schneider et al. 
1992) 



AGMD LS 



so that the lens equation is simply rewritten as 
where /3 = /3/0 E , 9 = O/0 E , and 



1 i,k=i 



Qij 



3=1 



0| 4 



(4) 



(•») 



(6) 



Here, the renormalised quadrupole moment is defined as 



Qi 



c 2 D s Q t , 



2GM*D L D LS 



(7) 



Eq. JSJl is a set of coupled equations which are seventh-order 
in 8 1 and 8 2 , respectively. It is apparent that we can hardly 
treat it analytically. 

Before examining Eq. JSJ in detail, we estimate the or- 
der of magnitude of the renormalised quadrupole moment. 
We can put \Qij\ ~ eMR 2 , where s denotes a distortion 
parameter and R is a typical size of the lens object. Trun- 
cations up to quadrupole moment in multipole expansions 
are valid if the asymmetry of the system is sufficiently small, 
say e < 1/10, though we don't perform any expansion in e. 
The renormalised quadrupole moment is rewritten as 



2s 



R 
Rv 



(8) 



where R-q = is the Einstein ring radius. In most 

cases, R is shorter than R E . Thus, \Qij\ must be less than 
2e, which is usually 1/5 at most. 

For instance, we assume a nearby solar-type star at 
lOpc, namely R E ~ 10 7 km, R ~ 10 6 km. If it rotates much 
faster than the Sun, say with the rotational surface velocity 
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on the equatorial plane v ~ lOkm/s, the oblateness due to 
the centrifugal flattering is of the order of 



10' 



V M A: 



J? 



10 6 km/ V lOkm/s 



(9) 



Then, Qij becomes about 10~ 5 . If a light ray passes near the 
surface of the star, the leading term of the deflection angle 
is about 1 arcsec. and the correction due to the quadrupolc 
moment can be as large as 10 microarcsec, which will be 
detectable in principle by a near future astrometry mission 
such as GAIA and SIM, though such a probability may be 
extremely small. 

Henceforth, we adopt a frame in which the reduced 
quadrupole moment is diagonalised as 



Qi 



e 
-e 



(10) 



Without loss of generality, we can assume e > 0. By denoting 
f3 — (a,b) and = {x,y) = (r cos <f>, r sin (f>) for r > 0, it is 
convenient to rewrite Eqs. J5J and JSJ into a matrix form 



where 



cos 4> cos 3(f> 
sin <j> sin 3<p 



f 



e 



(11) 

(12) 
(13) 



Here, / and g depend on the mass and quadrupole moment, 
respectively, where we should remember that our renormal- 
isation was done in units of the Einstein ring radius so that 
the mass does not appear explicitly. 



2.2 Sources on the symmetry axes 

First of all, we consider a source located at the origin (a, 6) = 
(0,0), namely exactly behind the lens object. For 6 = 0, 
the second line of Eq. 1 1 1 li means (1) sin<^ = or (2) / + 



#(3 -4 sin 2 



0. The case (1) means y = 0, consequently 



r — x . Then, the first line of Eq. 1111 becomes 
X 4, - x 2 - e = 0, (14) 
where we used r / 0. This equation is solved immediately 



as 



x = ± 



1 + Vl + 4e 



(15) 



where we used \fl + 4e > 1. Next, we consider the case (2), 
whose condition is rewritten as 



/ + 4gcos <j> = g. 



(16) 



By substituting this into the first line of Eq. (II IB . we obtain 
gcostj) = 0. Because of j / 0, we find cos</> = 0, which 
means x = and r 2 = y 2 . Then, Eq. (1161 becomes simply 
/ = g, which is rewritten as y 4 — y 2 + e = 0. This is solved 

as 



'1± VI ~4e 



(17) 



which are real if and only if e < 1/4. Hence, we find two 
images given by Eq. 1151 on a;-axis, and four images by Eq. 

on y-axis if e < 1/4, or otherwise only the first two 
images on a;-axis. 

Now, let us consider a case that a source is located 
on a symmetry axis but off the origin. For simplicity, we 
assume 6 = and a 7^ 0. The case of a — and 6/0 
is obtained simply by transforming e — > — e and x «-» y as 
discussed later. For 6 = 0, the second line of Eq. (I1H means 
(a) sin0 = or (b) f + g(3 - 4sin 2 <j>) = 0. The case (a) 
means y = and r 2 = x 2 . Then, the first line of Eq. JTT1 
becomes 



u(x) 



e = 0, 



(18) 



where we used r/0. The discriminant for this fourth-order 
polynomial (van der Waerden 1966) is computed as 



Da 



-e(16 + 128e + 256e 2 + 4a 2 + 144a 2 e + 27a 4 e), (19) 



which is negative. This implies that Eq. (I18H has only two 
real roots. Furthermore, we find tt(0) = — e < 0. The corre- 
sponding two images are located on x-axis, one positive and 
one negative each. 

The case of a = and 6 ^ is more complicated. The 
transformation mentioned above gives us 



v(y) = y -by -y + e = 0, 



(20) 



which is the second line of Eq. (II IB . For this equation, the 
discriminant becomes 



Z> 4 = e(16 - 128e + 256e 2 + 46 2 
D4 — has two roots 



where we defined 



1446 2 e - 276 4 e) 



(21) 
(22) 



6± = 3^[(1 ~ 2 2 3 2e ) ± (1 + 2 2 3e) 3/2 ]. (23) 

One can show that b 2 _ < and b\ > 0. Hence, for 6 2 > b\, 
we find D4 < 0, which implies two real roots, namely two 
on-axis images. On the other hand, for 6 2 < b\, we find 
D4 > 0, which implies either four or no real roots. Further 
investigations are thus needed. We use the derivative of v(y) 
as 

v'(y) =y(V- 362/- 2). (24) 
The solutions of v'(y) = are y = 0, y±, where we defined 

(25) 



y± = ^(36± V3 2 6 2 + 2^ 



If v(y+) < and v(y~) < 0, four real roots exist, while no 
real roots for v(y+) > and v(y-) > 0. In the case of D4 > 
0, the former condition is equivalent to v(y+) + v(y-) < 0, 
which is solved as 

(26) 



6 2 >6 2 = < L{-l + Z- 1 ' i y/T+¥'e) 



The term 6 2 is positive if and only if e > 1/4. For e > 1/4, 
one can show 6 2 > 6+. If e > 1/4, we thus find no images 
for 6 2 < b\ and two for for 6 2 > b\, while for e < 1/4, four 
images exist if 6 2 <b\, and two images if 6 2 > b\. 

We consider the case (b), whose condition is rewritten 

as 
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/ + 45COS <f> = g. 



(27) 



By substituting this into the first line of Eq. 1111 . we obtain 
a = —2gcos(j). Because of g 7^ 0, we find 



ar 
~~2e 



(28) 



Substitution of this into Eq. Ij27^ gives us a cubic equation 
for r 2 as 



2 6 4.2 2 n 
a r — er + er — e = 0. 



(29) 



Because of the symmetr y with resp ect to x-axis, images may 
appear at (r cos <f>, ±r^/l — cos </> 2 ), where r is a solution of 
Eq. and cos</> is obtained by substituting the solution 
for r into Eq. 

Here, we should note a constraint on r, which comes 
from Eq. @) for | cos 0| < 1. The constraint is expressed as 



r < 



1/3 



(30) 



One can investigate the number of off-axis images ac- 
cording to Eqs. 129H and (1301) . The detailed discussion is 
given in the Appendix A. It shows that the maximum 
number of off-axis images is four if < e < 1/4, six if 
1/4 < e < 1/3, and two if e > 1/3. Similar results in a case 
of a — and b 7^ are given in the Appendix B. 



2.3 Off-axis sources 

Here, let us consider off-axis sources, n/0 and b 7^ 0. Then, 
Eq. 1111 implies sin cj> 7^ and cos 7^ 0, which are equivalent 
to sin 2<f) 7^ 0. By using an inverse matrix, which always 
exists if sin 2<j> 7^ 0, Eq. [11 1 1 is rewritten as 



sin 2(f> 



sin 2>4> 
— sin d> 



- cos 2>4> 
cos d> 



(31) 



Thus r-dependent parts are separated from ^-dependent 
ones. A key is the first line of Eq. 1311 . which is re-expressed 
as r 2 = Fr + 1. By using this recursively, we obtain r 3 = 
(F 2 + l)r + F. Substituting this into r 3 in the left hand side 
of the second line of Eq. 1311 . we obtain 



FG - 



G(F 2 + 1) 



(32) 



whose right-hand side depends only on <f). This equation im- 
plies that G does not vanish for e/fl because r is finite. We 
should note that r cannot be expressed as a function only 
of cot 4> but also of cos cf>. 

Now, we rewrite F and G as 



F = 

G = 

where 

u = 
F = 



F 



2(1 +v?) cose 
bu — a 



2 cos ( 



cot < 



a(3u — 1) — bu(u — 3). 



(33) 
(34) 



(35) 
(36) 



Thus, by using Eqs. (13211 . I33H and 134H . we can express (a;, y) 
as a function only of u, 



F 



2(1 



F< 



+ 



2eu 2 



V = 



4u 2 (l + u 2 ) 
Fp 

2u(l+u 2 ) 



(bu-a)(l 
4- J ' 



w 2 ) 



+ 



2eu 



FS 



4u 2 (l + u 2 



(bu-a)(l + u 2 ) 
^ 1 
+ 1 X 



(37) 



(38) 



This fact implies that there must exist a master equation 
determining u. Substituting Eq. 1)32^ into r in the first line 
of Eq. I3H , we obtain a tenth-order equation for it = cot 4> 



10 

Cku k 

fc=0 



0. 



(39) 



where each coefficient is a polynomial in a, b and e as 



C10 
eg 

C8 
C7 
C6 

C5 
C4 

C3 
C2 
Cl 
Co 



, 4 

e, 



-10a6 3 e, 



46 2 + 36 2 e(4 + 12a 2 - 
-8a6-6a6e(8 + 9a 2 



■36"), 
- lib 2 



(40) 
(41) 
(42) 
(43) 



4a 2 +i 



b 2 +3e(12a 2 + 9a 4 



8b 2 



52aV + 96 4 ) 



-16e 2 , 

-16a6+ 126abe(a 2 - 6 2 ), 

8a 2 + 46 2 + 3e(8a 2 - 9a 4 - 126 2 

-16e 2 , 

-8ab + 6abe{8 - 11a 2 + 96 2 ), 



= 4a" 



3a 2 e(4 



3a 2 + 12fo 2 ), 



10a 3 &e, 



(44) 
(45) 
9b 4 ) 
(46) 
(47) 
(48) 
(49) 
(50) 



As seen from Eq. 1391 , u = cot <f) neither vanishes nor 
diverges because C10 7^ and Co 7^ for a, b and e/0. The 
number of images is constant unless a source crosses the 
caustics. Therefore, the maximum number of real solutions 
for Eq. 1391 is eight, which equals that for an on-axis source 
as shown in the preceding subsection. One can see also that 
there is a <-+ cio-i correspondence with a <-> 6 and e <-> 
— e transformation. This is taken as a justification for the 
statement regarding the relation between (a = 0, b 7^ 0) and 
(a/ 0,6 = 0). 



3 DISCRIMINANT AND CAUSTICS 

As well-known in algebra, a discriminant for a real polyno- 
mial tells us where the number of real roots changes (for 
instance, van der Waerden 1966). A discriminant for Eq. 
1391 is obtained for instance by using Mathematica (Wol- 
fram 2000) as 

D10 = 2 42 a 12 b 12 e 6 W 2 K, (51) 



where 
W = 



(a 2 + 6 2 ) 2 + (a 2 + 6 2 ) 3 -2e(a 2 
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+ e 2 [l-(a 2 + 6 2 )+3(a 2 + & 2 ) 2 ] 
-e 4 [2-3(a 2 + & 2 )]+e 6 , (52) 
and K is a lengthy polynomial of 494 terms as 



K = 2 2 3 lb a V(2 + 18e + 276 e) + 



(53) 



This discriminant Z?io gives us a condition for changes in the 
number of real roots for Eq. 1391 . corresponding to changes 
in the number of images. 

It should be noted that W never vanishes as shown 
below, which means that r is uniquely determined by Eq. 
13211 . This is in contrast to a case of a binary gravitational 
lens, for which a squared term such as W 2 can vanish so that 
a function of r can take a form of "0/0" for particular cases 
(Asada et al. 2002). By paying attention to an asymmetric 
term proportional to a 2 — b 2 in Eq. 1521 . we rewrite W as 

W = (a 2 + b 2 - e + e 2 ) 2 [a 2 + b 2 + (1 + e) 2 ] + 46 2 e, (54) 

which shows that W is always positive for b ^ and e > 0. 
The caustics are thus expressed simply as K — 0. 



4 CRITICAL CURVES 

In the preceding section, the explicit expression for the caus- 
tics has been obtained. If we map the expression onto the 
lens plane by the lens equation, we can obtain the corre- 
sponding expression for the critical curves. In practice, how- 
ever, it is a quite difficult task. Hence, let us study the crit- 
ical curves in a different way. By a straightforward calcula- 
tion of the vanishing Jacobian of the lens mapping by Eq. 

CD, 

d(a,b) _ 1 d(a, b) 
d(x,y) rd(r,4>) 

we obtain a simple expression for the critical curves as 
r 8 = r 4 + 9e 2 + 6er 2 cos 2(j>. 



0. 



(55) 



(56) 

This curve is classified topologically into two as follows. Be- 
cause of | cos2<^)| < 1, Eq. 1561 gives a constraint on r 

(57) 
(58) 



which is equivalent to 
r 4 < jr 2 + 3ej, 



and 



r > \r 



3e|, 



(59) 



where we used r 4 > 0. By solving carefully Eqs. 1581 and 
EH, we find for e > 1/12 



-l + VT+Ui ^ 2 / l + 71 + 126 
< r < , 



and for e < 1/12, 



-1 + VT+m a 
< r < 



1 - VI - 12e 



1 + VI - 12e ^ 2 ^ 1 + VI + 12e 
- < r < - 



(60) 



(61) 



(62) 



2 - - 2 

This implies a topological feature of the critical curves: The 
curves are simply connected if e > 1/12, or otherwise they 
are not. 




Figure 1. Caustics for a lens with quadrupole moment e = 0.05. 
Sources are located at (0, 0), (0.05, 0) and (0.2, 0), denoted by the 
triangle, filled disk and square, respectively. 



When Eq. 1)56^ is treated as a ^-parameter representa- 
tion of the critical curves, we must solve it as a fourth-order 
equation for r 2 , so that r 2 can be multi- valued. As a conse- 
quence of r > 0, furthermore, some domains of the parame- 
ter (j> & [0, 2-7r) can be excluded. It seems thus much simpler 
to consider Eq. (1561 as a representation with a parameter r 
which is allowed by Eq. (1601 for e > 1/12 or by Eqs. 1611 
and 1621 for e < 1/12. The representation becomes 



(x,y) = (±7-4/^4— ,±tm "' '' 



where we defined 
r 4 - 9e 2 



h = 



6er 2 



(63) 



(64) 



In addition, substitution of this parametric representa- 
tion into the right-hand side of Eg. (1111 gives us a represen- 
tation of the caustics with the same parameter as 



a(r) 



b(r) 



±[f + g(2h-l)] 



±[f + g(2h + l)} 



1 + h 



1 - h 



(65) 



(66) 



As illustrations, the caustics and critical curves for 
e = 0.05 < 1/12 are given by Figs. 1 and 2, respectively. 
Actually, there are three closed loops of critical curves. A 
case of e = 0.1 > 1/12 is shown by Figs. 3 and 4. 



5 CONCLUSION 

We have re-examined the lens equation for a gravitational 
lens due to a deformed star modeled by using quadrupole 
moment. First, we reduce the lens equation to a single real 
tenth-order polynomial. Consequently, an analytic expres- 
sion for the caustics is given by Eq. I|53fl as K — 0, though 
it is too lengthy to handle. The critical curves are simply 
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Figure 2. Critical curves for a lens with quadrupole moment e = 
0.05. The images correspond to the sources in Fig. 1. Triangles, 
Filled disks and squares are overlapped around (0, ±0.2). 



Figure 4. Critical curves for a lens with quadrupole moment 
e = 0.1. The images correspond to the sources in Fig. 3. Triangles, 
Filled disks and squares are overlapped around (0, ±0.3). 



1 




-1 [ J 

-1 -0.5 0.5 1 

a 

Figure 3. Caustics for a lens with quadrupole moment e = 0.1. 
Sources are located at (0,0), (0.1,0) and (0.3,0), denoted by the 
triangle, filled disk and square, respectively. 



expressed as Eq. 11561 or Eq. 16311 . The parametric represen- 
tation of the caustics is given by Eqs. 1651 and 1661 . Hence, 
the present result must be helpful for understanding the 
gravitational lensing due to a deformed star. In particular, 
it gives a nice approximation to a lensing study based on a 
realistic numerical simulation of a compact object. 

For instance, the image position is unstable for a source 
close to the caustics, so that careful numerical computations 
are needed. The present approach will make it easy to study 
such a case. This is a future subject. 
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APPENDIX A: CLASSIFICATIONS OF 
OFF-AXIS IMAGES I. 

Let us investigate the number of off-axis images according 
to Eqs. ((23 and 

in a case of a 7^ and 6 = 0. They are 

rewritten as 



g(p) 

P< 



2 3 
a p 



2/3 



ep + ep — e 



: Pa, 



0. 



where we defined p — r 2 . The discriminant of Eq. IA1I 



2-3 2 e)a 2 -(e-2 2 e 2 )]. 



D s = -e 3 [3 3 ea 4 + (2 
D3 = has two roots as 
a 2 = A±, 
where we defined 

A± = ^[(3 2 e-2)±2(l-3e) 3/2 ] 

We find A_ > if e > 1/4, and otherwise A_ < 0. We find 
g(po) as 



(Al) 
(A2) 

I is 

(A3) 

(A4) 
(A5) 



9(Po) = £ ( 3el/3fl4/3 



2 4/3 e 2/3 + 2 2/3 a 2/3N 



The condition that g(po) > is 
2 



a > A 



3 :i f: 



(1 + 3 2 e) + (1 + 3e)Vl + 2 2 3e], 



(A6) 



(A7) 



where we used the positivity of a 2 . 

In order to analyse the number of the valid solutions 
satisfying Eq. 1A2L one needs to investigate the signature 
not only of g(po) but also of g'(po). In a case of g(po) < 
and g'(po) > 0, for instance, one can show that there exists 
a root p > po, which corresponds to an unphysical solution 
r violating the condition (1301 . The condition that g'(po) < 
is 



a < Ai 



2 5 3 3 e 



[-( 1 + 2 3 3 2 e) + ( 1 + 2 3 3e) y/l + 2 5 3e] , (A8) 



where we used the positivity of a . 

The number of the valid solutions are found by figur- 
ing out the relative order among a 2 , A±, Ao and Ai, which 
depends on the value of e. Straightforward but tedious com- 
putations give the following result. 
(1) 0<e< \ 

(la) A_ < < a 2 < A < A + < A x 

There are three positive roots. Only two of them are valid 
and thus four off-axis images exist, 
(lb) A_ < A < a 2 < A + < A 1 

There are three positive roots. Only one of them is valid and 
thus two off-axis images exist, 
(lc) A_ < A < A+ < a 2 < Ai 

There is only one positive root, which is still valid. Thus two 

off- axis images exist. 

(Id) A- < Ao < A+ < Ai < a 2 

There is only one positive root, which is still valid. Thus two 
off-axis images exist. 



(2) \ < e < £ 

The relative order between A- and Ai is not unique in this 
case. 

(2a) < a 2 < A- , Ai < Ao < A+ 

There is only one positive root, which is not valid. Thus no 

off-axis images exist. 

(2b) < A_ < a 2 < Ai < A,, < A+ 

There are three positive roots. Only two of them are valid 
and thus four off-axis images exist. 
(2c) < Ai < a 2 < A- < Ao < A+ 

There is only one root, which is not valid. Thus no off-axis 
images exist. 

(2d) < A_, Ai < a 2 < A < A+ 

There are three positive roots. Only two of them are valid 
and thus four off-axis images exist. 
(2e) < A_, Ai < A < a 2 < A+ 

There are three positive roots, all of which are still valid. 
Thus six off-axis images exist. 
(2f) < A_, Ai < A < A+ < a 2 

There are only one positive root, which is still valid. Thus 
two off-axis images exist. 

(3) A < c < 1 

(3a) < a 2 < Ai < A < A_ < A + 

There is only one positive root, which is not valid. Thus no 

off- axis images exist. 

(3b) < Ai < a 2 < A < A_ < A+ 

There is only one positive root, which is not valid. Thus no 

off- axis images exist. 

(3c) < Ai < A < a 2 < A_ < A+ 

There is only one positive root, which is still valid. Thus two 

off-axis images exist. 

(3d) < Ai < A < A_ < a 2 < A+ 

There are three positive roots, all of which are still valid. 
Thus six off-axis images exist. 
(3e) < Ai < A < A_ < A+ < a 2 

There is only one positive root, which is still valid. Thus two 
off-axis images exist. 

(4) § < e 

(4a) < a 2 < Ai < A 

There is only one positive root, which is not valid. Thus no 
off-axis images exist. 
(4b) < Ai < a 2 < A 

There is only one positive root, which is not valid. Thus no 
off- axis images exist. 
(4c) < Ai < A < a 2 

There is only one positive root, which is still valid. Thus two 
off-axis images exist. 



APPENDIX B: CLASSIFICATIONS OF 
OFF-AXIS IMAGES II. 

We classify the number of off-axis images in a case of a = 
and b 7^ 0. The transformation of a <-> b and e <-> — e changes 
Eq. (J53J as 

6 2 r 6 + er 4 - er 2 - e 2 = 0. (Bl) 
The constraint by Eq. 1301 is also transformed into 



r < 



1/3 



(B2) 
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Let us investigate the number of off-axis images accord- 
ing to Eqs. JBTJ and JH1J in a case of a = and 6 7^ 0. They 
are rewritten as 

%) = b 2 q 3 + eq 2 -eq-e 2 = 0, (B3) 



q< 



2/3 

= qo, (B4) 



where we defined q = r . 

One can show that the positive root is always only one 
because h(0) < and h'(0) < 0. This situation is much 
simpler than that of a 7^ and b = 0. Namely, one can show 
that the valid solution is still one if h(qo) > 0, while we have 
no valid solution if h(qo) < 0. We find h(qo) > always 
holds if e > 1/12. If e < 1/12, h(q ) > is equivalent to 
b 2 < B- or b 2 > B+, where we defined 

2 



B± = —[(1 - 3 2 e) ± (1 - 3e)Vi - 2 2 ■ 3c]. (B5) 

As a result, the number of the off-axis images is classified as 
follows. 

(1) e<i 

(la) < b 2 < B- 

There is only one positive root, which is still valid. Thus two 
off-axis images exist, 
(lb) B- < b 2 < B+ 

There is only one positive root, which is not valid. Thus no 
off-axis images exist, 
(lc) B+ < b 2 

There is only one positive root, which is still valid. Thus two 
off-axis images exist. 

(2) e > i 
(2a) < b 2 

There is only one positive root, which is still valid. Thus two 
off-axis images exist. 
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